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ABSTRACT 

We explore the effects of the deflagration to detonation transition (DDT) density on the production of ^^Ni 
in thermonuclear supernova explosions (type la supernovae). Within the DDT paradigm, the transition density 
sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear 
statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled sta- 
tistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions 
that can, with a particular choice of transition density, produce a similar average and range of ^''Ni masses 
to those inferred from observations. Within this framework, we utilize a more realistic "simmered" white 
dwarf progenitor model with a flame model and energetics scheme to calculate the amount of ^^Ni and NSE 
material synthesized for a suite of simulated explosions in which the transition density is varied in the range 
1-3 X 10^ g cm"^. We find a quadratic dependence of the NSE yield on the log of the transition density, which 
is determined by the competition between plume rise and stellar expansion. By considering the effect of metal- 
licity on the transition density, we find the NSE yield decreases by 0.055 ± 0.004 Mq for a 1 Zq increase 
in metallicity evaluated about solar metallicity. For the same change in metallicity, this result translates to a 
0.067 ± 0.004 Mq decrease in the ^''Ni yield, slightly stronger than that due to the variation in electron fraction 
from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat 
ambiguous, but the dependence we find is comparable to that inferred from some studies. 

Subject headings: hydrodynamics — nuclear reactions, nucleosynthesis, abundances — supernovae: general 
— white dwarfs 



1. INTRODUCTION 

Type la supernovae are bright stellar explosions that are 
characterized by strong P Cyg ni features in Si and by a lack of 
hydrogen in their spectra (see Hillebrandt & Niemeveij|2000t 
[pilippenko 1997 and references therein). Obser vations of type 
la supernovae (se r ving as distance indic ators (lPhinipslll993l: 
iRiess et all 119961: lAlbrecht et all l2006h ) are at present the 
most po werful and best proved technique for studying dark 
energy (Riess et al."1998; 'Perlmutter et al."1999!; iKolb et all 
12006; Hicken et al. 2009; Lampeitl et al. 2010), and, accord- 
ingly, there are many observational campaigns underway 
striving to gather information about the systematics of these 
events an d to better measu re the expansion history of the Uni- 
verse (see ' Kirshnerll2009l and references therein). 

The most widely accepted model for these events is the 
"single-degenerate" scenario, which is the thermonuclear ex- 
plosion of a white dwarf composed principally of '^C and 
'^O that has accreted m ass from a companion (for a review 
of explosion models, see ' Hillebrandt & Niemeverl2000l;lLiviol 
HOOO; Ropke 2006). The peak brightness of the supernova is 
set by the radioactive decay of ^^Ni produced in the explosion 
to ^^Co and then to ^''Fe. The empirical "brighter is broader" 
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relation between the peak brightness of the light curve and 
the decay time from its maximum is understood to follow 
from the fact that both the luminosity and opacit y are set by 
the mass of ^^Ni svnthesi zed in the explosion (,ArnetJll982l; 
iPinto & Eastman 2000; K asen & Wooslevll2007h . Because of 
the dependence of the light curve on the amount of radioac- 
tive ^^Ni synthesized during an explosi on and the ab ility to 
infer the ^^Ni yield from observations (JMazzali et al.l l2007). 
research into modeling thermonuclear supernovae typically 
focuses on the production and distribution of ^^Ni as well as 
other nuclides (such as -**Si) as the measure with which to 
compare models to observations. 

One-dimensional simulations of the single-degenerate case 
showed that the most successful scenario is an initial de- 
flagration (subsonic reaction front), born in the core of the 
white dwarf (WD) which at some point becomes a (su- 
personic) d etonation, i.e. a deflagration-detonation tr ansi- 
tion (DDT, Kh okhlovl [19911 iHoflich & KhokhlovllT996h . A 
delayed det onation naturally acco unts for the high-velocity 
Ca features jKasen & Plewall2005b and the chemical stratifi- 
cation of the ejecta. While these ID models are able to repro- 
duce observed features of the light curve and spectra, much 
of the physics is missing from these models. The presence of 
fluid instabilities during the deflagration warranted the devel- 
opment of multidimensional models, allowing a physically- 
motivated calculation of the velocity of the burning front and 
thus removing a free parameter. By relaxing the symmetry 
constraints on the model, buoyancy instabilities are naturally 
captured leading to a strong dependence on the initial condi- 
tions of the deflagration that often result in too little expan- 
sion of the star by the ti me DDT conditions used in previ- 
ous IP stud i es are met ( Nierneveret~anil996t ICalder et^ 
I2003L 120041: iLivne et al.1 [20051). Mult"idimensional models 



may reach the expected amount of expansion prior to the 
DDT with the choice of particular ignition conditions and 
thus retain the desirable features from one-dimensional mod- 
els (Gamezo et al. 2004; Plewa et al. 2004; Ropke et al. 2006; 
[Jordan et al . 2008). 

Due to the strong dependence on ignition conditions, multi- 
dimensional simulations of the DDT model are able to pro- 
duce a wide range peak luminosities (via the production of 
a range of ^^Ni yields) consistent with a common exp losion 
mode suggested from observations (iMazzali et al.l2 007). Dif- 
ferences in the mass of synthesized ^^Ni can follow from 
properties such as metallicity and central density of the pro- 
genitor and/or differences in the details of the explosion 
mechanism such as the density at which the transitio n from 
deflagration to detonation occurs. Timm es et al.l (l2003h found 
that metallicity should affect the ^^Ni yield based on ap- 
proximate lepton number conservation. The metallicity sets 
the fractional amount of material in nuc lear statistical equi- 
librium (NSE) that is radioactive ^^Ni. iBravo et al] (1201 Ol) 
found a stronger dependence on metallicity due to a signif- 
icant amount of ^^Ni that is synthesized during incomplete 
Si -burning. 

Observational results to date are consistent with a shal- 
low dependence of ^^Ni mass on metalhcity but are un- 
able to rule out a ti-end entirely (Gallag her et al.ll2005L 120081; 
[Neill et al. 2009; Howell et al. 2009). Determining the metal- 
licity dependence is challenging not only because the effect 
appears to be small, but also due to the difficulty in measur- 
ing accurate metallicities for the parent stellar population and 
problems with strong systematic effects associated with the 
mass-metallicity relationship within galaxies (Gallagher et al. 
|2008|) . This effect is also difficult to decouple from the ap- 
parently stronger effect of the age of the parent stellar popula- 
tion on the rnean brightness of SNe la ( [Gallagher et al. 2008; 
iHowell et al.ll2009l;lKrueger et alj |2010^. Ho well et alj (i2009i) 
note that the scatter in brightness of this observed relation is 
unlikely to be explained by the effect of metallicity. In gen- 
eral, the source of scatter can be explained by the develop- 
ment of fluid instabilities during the deflagration phase that 
contribute to differing rates of expansion between instances 
of supernovae. However, by considering the effect of metal- 
licity on the DDT density, scatter in the metallicity relation 
may be enhanced beyond its intrinsic value inferred from fluid 
in stabilities. 

iTownslev et al.l (l2009l) investigated the direct effect of 
metalhcity via initial neutron excess and found it to have a 
negligible influence on the amount of material synthesized to 
NSE. However, the neutron excess sets the amount of material 
in NSE that favors stable Fe-group elements over radioactive 
^^Ni. Therefore, the initial metallicity directly influences the 
yield of ^*Ni. In this work, we expand that study to include 
the potential indirect effect of metallicity in the form of the 
^^Ne mass fraction (X22) through its influence on the density 
at which the DDT takes place. To this end, we explore the ef- 
fect of varying the transition density, a proxy for varying the 
microphysics that determine the conditions for a DDT. The 
conditions und er which a DDT occu r s are still a subject of 
debate . Niemever & WooslevI (Il997h . iNiemeyer & KersteinI 
(I1997I) . and Khokhl ov et al.l (Il997h proposed that a necessary 
condition is the transition to a distributed burning regime, 
in which turbulence disrupts t he reaction zone of the flame. 
More recent numerical studies (iPan et al. "2008; Woosley et al. 
12009; Schmidt et al. 2010 ) that describe the conditions for a 
DDT include dependencies on the turbulent cascade and the 



growth of a critical mass of fuel with sufficiently strong tur- 
bulence. For this study, we assume the DDT density to be 
the density at which thermonuclear burning is expected to en- 
ter the distributed regime. This choice links the explosion 
outcome to the dynamical evolution of the progenitor density 
structure during the deflagration phase. By analyzing the ef- 
fect of transition density on the NSE yield, we can later ana- 
lyze how the details of the microphysics affect the DDT den- 
sity. For the purposes of this study, we will consider only the 
effect of X22 on the DDT density. In reality, the '^C mass 
fraction will also be important in determining the DDT den- 
sity, but we choose to leave the exploration of the effect of '^C 
to future work. Many other possible systematic effects exist 
that are outlined in Townslev et al. (2009), such as the central 
ignition density dKrueger et alJ2010 ), which are all held fixed 
in this study. The interdependence of all of these effects will 
be considered in the construction of the full theoretical picture 
in a future study. 

We describe the details of our model in ^and ^ and the 
properties of our statistical sample in ^ We present our find- 
ings on the dependence of transition density on NSE yield 
in ^5.11 In ^5.21 and ^5.31 we assume a dependence of the 
transition density on ^^Ne content and construct the functional 
dependence of the ^^Ni yield on metallicity through the ^^Ne 
content. In ^ we discuss our conclusions and future work. 

2. PARAMETERIZED REALISTIC WHITE DWARF PROGENITOR 

In order to include relevant processes in explosion mod- 
els, we first estimate the compositional profile of the progen- 
itor white dwarf just prior to the birth of the flame. We be- 
gin by estimating the compositional profile resulting from the 
evolution of the post-main-sequence star that later becomes 
a white dwarf. Recall that the initial metallicity of the star 
is, by mass, almost entirely in the form of CNO. As a result 
of the C-N-O cycle, these all end up in helium layers as '"^N, 
the target of the slowest proton-capture in the cycle. Subse- 
quent reactions during helium burning convert '"*N into ^^Ne; 
ther efore, X2?, is proport ional to the initial metallicity of the 
star dTimmes et al.ll20()3l) . The composition of the inner por- 
tion of the star (sa 0.3-0.4 Mq) is set during core helium 
burning resulting in a reduced carbon mass fraction with re- 
spect to that of the outer layers, wh ich is set by shell bur ning 
on the asymptotic giant branch (see lStraniero et all2003L and 
references therein). 

At some point after the white dwarf has formed, it begins 
to accrete material from its companion. As the mass of the 
accreting white dwarf approaches the Chandrasekhar mass 
limit, the core temperature and density increase such that car- 
bon begins to fuse. The energy released by carbon burning 
drives convection in the core. The convective carbon burn- 
ing ("simmering") phase lasts approximately 1000 years be- 
fore the central t emperature is high en ough to spark a ther- 
monuclear flame dWooslev et alj|2004h . During the simmer- 
ing phase, carbon is consumed from the convective core. Con- 
currently, though, the convective zone grows with increasing 
central temperature, pulling in relatively carbon-rich material 
from the outer layer. Thus, the net effect is to increase the car- 
bon abundance in the convective region. We show the growth 
of the convective zone in Figure [H where dashed lines show 
the compositional profile prior to the simmering phase and 
the corresponding solid lines show the compositional profile 
at the start of our simulations of the explosion. 

For our WD models, we consider a parameterized 3-species 
compositional structure consisting of ^^C, ""O, and ^^Ne, 
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TABLE 1 

Composition of the core and outer layer of the progenitor 

white dwarf. 









Core 




Outer Layer 




Xn 


^16 


xu 


Mc (Mq) 


Xn 


X[6 


X,', 


Pre-Simmering 
Pre-Deflagration 


0.22 
0.40 


0.76 
0.57 


0.02 
0.03 


0.30 
1.16 


0.50 
0.50 


0.48 
0.48 


0.02 
0.02 



which is sufficient to capture the nuclear burning rates. Since 
^^Ne is the only element in our parameterization that has more 
neutrons than protons, the neutron excess from simmering is 
accounted for by a parameterized ^^Ne mass fraction 



X^2=^22 + A^22(Ay,), 



(1) 



where X22 is proporti onal to the initial m etallicity coming 
from hehum burning dTimmes et alj 120031) and AX22(AYe) 
represents the effective enhancement of ^^Ne in the core fol- 
lowing from the change in electron fraction during convective 
carbon burning. The electron fraction (Yg) is related to X22 by 



22^22+2(1-^22) 



(2) 



such that a change in X22 constitutes a change in Yg. The prime 
on X20 in Equation ^ indicates inclusion of effects of neu- 
tronization during carbon burning and this convention applies 
to expressions below. We choose the composition of the WD 
at the end of the simmering phase to consist of X12 = 0.4 and 
X22 = 0.03 in the core and X12 = 0.5 and X22 = 0.02 in the outer 
layer. Note that X22 = ^22 ^'^ ^^^ outer layer because neutron- 
ization due to carbon burning only occurs in the convective 
core. 

For comparison, a compositional profile can be estimated 
for the WD prior to the onset of carbon burning. For sim- 
plicity, consider the production of one neutron for every 6 '^C 
burned (i.e., dYe/dYi2 ~ 1/6, where F12 is the molar abun- 
dance of '^O (IChamulak et al.ll2008t IPiro & Bildstenll2O08h . 
For AX22 = 0.01 in our progenitor model, this constitutes 
burning 0.04 Mq of carbon during simmering. Assuming that 
prior to simmering, the core mass is w 0.3 Mq dStraniero et alJ 
l2003h and the outer layer consists of X12 = 0.5 material, we 
can account for all '^C and conserve the total mass of the WD 
to estimate X12 ~ 0.2 in the carbon-reduced core of the WD 
prior to simmering as shown in Figure [T] 

Table [T] shows the composition of the progenitor white 
dwarf before and after the simmering phase using the param- 
eterized ^^Ne mass fraction as well as the core mass. Note 
that throughout this study, we choose to neglect any variation 
of the amount of neutronization during simmering due to the 
initial Z. Accordingly, AX22(AYe) is treated as a constant and, 
therefore, dX'22 = dX'22 and derivatives involving the true ^^Ne 
mass fraction proportional to metallicity, X22, are equivalent 
to derivatives involving the parameterized ^^Ne mass fraction, 

^22- 

Because we c onsider the enhanceme nt of the laminar flame 
speed by ^^Ne (IChamulak et al.ll2007h . we need to consider 
the effects of our parametrization of the ^^Ne content with 
care. In actuality, the neutronization during carbo n simmering 
produ ces '^C, ^^Na, and ^"Ne and not -^Ne (Cha mulaket alJ 
l2008h . A priori, this difference could alter the nuclear burn- 
ing rates and hence the laminar flame speed and width. Ta- 
ble|2]shows the laminar flame speeds and widths for the same 



TABLE 2 
Flame speeds and widths changing -^Ne to carbon-simmering 

ASHES holding Ye = 0.498636 FIXED AND USING X12 = 0.4. PRIMED 

QUANTITIES PARAMETERIZE THE EFFECTS OF NEUTRONIZATION 

DURING CARBON BURNING AS ADDITIONAL ^^Ne. 



P9 (g cm-3) 


.s'(kms-') 


s (km s ' ) 


Diff. (%) 


0.1 


0.926 


1.012 


8.5 


0.2 


4.194 


4.570 


8.2 


0.3 


11.372 


12.106 


6.1 


0.4 


18.785 


19.466 


3.5 


0.5 


24.352 


25.057 


2.8 


0.6 


29.162 


29.916 


2.5 


0.7 


33.527 


34.322 


2.3 


0.8 


37.571 


38.401 


2.2 


0.9 


41.364 


42.228 


2.0 


1.0 


44.978 


45.871 


1.9 


P9(gcm-3) 


<5' (cm) 


5 (cm) 


Diff. (%) 


0.1 


1.4018 X 10^ 


1.3829 X 10^ 


-1.4 


0.2 


2.7439 X 10^ 


2.4097 X 10^ 


-13.9 


0.3 


8.5650 X 10* 


8.1055x10* 


-5.7 


0.4 


4.6600 X 10" 


4.9033 X 10* 


-5.0 


0.5 


3.5373 X Iff* 


3.4319 X 10* 


-3.1 


0.6 


2.6457 X 10" 


2.5756 X 10* 


-2.7 


0.7 


2.0980 X 10'* 


2.0437 X 10* 


-2.7 


0.8 


1.7167x10'' 


1.6820x10* 


-2.1 


0.9 


1.4557 X lO'* 


1.4199 X 10* 


-2.5 


1.0 


1.2562 X 10* 


1.2251 X 10* 


-2.5 



Ye exchanging AX22 for carbon-simmering ash for a compo- 
sition of 40% '^C, 2.0% ^^Ne, 55.5% '^O, 0.6% '^C, 0.9% 
2°Ne, and 1 .0% ^^Na by ma s s usin g the same method as de- 
scribed by iChamulak et al.l (120071) . The flame speeds and 
widths calculated using the parameterized ^^Ne mass fraction 
(3% by mass) are denoted with primes. For high densities 
(> 2.5 X 10^ g cm"-') in which the flame speed is not domi- 
nated by the buoyancy-driven Rayleigh-Taylor instability, the 
difference is < 5%; therefore, our parametrization of the neu- 
tronization via X22 accurately captures the corresponding en- 
hancement of the laminar flame speed in this regime. 

While the central temperature and central density of the 
progenitor just prior to the birth of the flame are primarily set 
by the accretion history of the white dwarf, which varies and 
is largely unknown, we choose a fiducial central density of 
Pc = 2.2 X 10^ g cm~^ and central temperature of r^ = 7 x 10^ 
K. We construct isentropic profiles of density and temperature 
in the (convective) core and isothermal profiles in the (ther- 
mally conductive) outer layer while maintaining hydrostatic 
equilibrium. We choose a fiducial isothermal temperature of 
the outer layer to be Tiso = 10^ K. The total mass of the WD 
progenitor model is 1 .37 Mq. 

Due to the difference in composition between the core and 
outer layer, requiring a neutral buoyant condition at this inter- 
face no longer reduces to a continuous temper ature as it does 
for a compositionally uniform white dwarf (iPiro & ChangI 
I2OO8I) . In reality, we might expect some convective overshoot 
and mixing between the core and the outer layer, but because 
the physics of convective overshoot are complex and not well 
understood, the width of the transition region is unknown. 
For simplicity, we assume no mixing region. The composi- 
tion, density, and thermal profiles of the progenitor used for 
this study immediately prior to deflagration are the solid lines 
shown in Figure [T] 

3. SIMULATION METHODS 




0.6 U.8 1.0 

Mass (M. ) 

Fig. 1 . — Composition, density and thermal profiles of the progenitor white 
dwarf star used in the simulations for this study (solid lines). The composi- 
tional profile of the progenitor prior to the simmering phase is also shown 
(dashed lines). 

We use a customized version of the FLASH E ulerian com- 
pressible adaptive-mesh hydrodynamics code tFrvxell et al.l 
l2000tlCalder et al.ll2002l) . Modifications to the code include a 
nuclear burning model, composition-dependent laminar flame 
speeds, particular mesh refinement criteria, and instructions 
for determining the conditions for a DDT (fo r details, see 
ICalderetal.ll2007t iTownsley et all l2007l 12009"). In particu- 
lar, the adaptive mesh refinement (AMR) capability of the 
code was utilized to achieve particular resolutions for burn- 
ing fronts (4 km) and the initial h ydrostatic star (16 km ) for 
which the solution has converged dTownslev et al.ll2009l) . 

The burning model is used for both subsonic (deflagra- 
tion) and supersonic (detonation) burning fronts. The lam- 
inar flame width for densities and compositions characteris- 
tic of a massive C-0 WD is unresolved on the scale of our 
grid (4 km, Chamulak et al. 2007); therefore, we use an artifi- 
cially thickened fla me represented by the advection-reaction - 
diffusion equation (lKhokhlovlll995l;IVladimirova et al.ll2006l) 
to which our nuclear energetics scheme is coupled. Appro- 
priate measures have been taken to ensure this coupling is 
acoustically quiet and stable such that the b uoyant instabil- 
ity of the burning front is accurately captured (iTownslev et al.l 
l2007h . The shock-capturing capabilities of FLASH naturally 
handle the propagation of detonation fronts given an accurate 
nuclear energetics scheme (Meakin et al. 2009). We note that 
the on l y co mmon components between our code and that of 
iPlewal (12007 ) are those publicly available as components of 
FLASH, which excludes all compo nents treating the nucle ar 
burning; differences are discussed in ITownslev et al.l (l2007h . 

In this section we discuss an improvement to our burning 
model and an improvement to the method by which a transi- 
tion from deflagration to detonation is made. The changes 
to the burning model reflect recent work to better match 
steady-state detonation structures that are partially resolved 
on the grid. This is important for obtaining accurate parti- 
cle tracks for post-processing nucleosynthetic yields. Addi- 
tionally, changes to the DDT method were necessary to en- 
sure consistency between individual simulations using differ- 
ent DDT densities and/or different initial configurations of the 
flame. 

3.1. Improved Burning Model 

For the calculations presented here we utilize the lat- 
est revision of a parameterized three-stage model for the 
nuclear burn i ng (ICalder et al.l 12007); ITownslev et all 120071; 
iMeakin et all l2009t ISeitenzahl et all l2009bt ITownslev et al.l 



|2009|) . The details of this latest version will be published 
separately (Townsley et al., in prep.) along with exten- 
sive comparisons to nuclear network calculations of steady- 
state detonations, but we summarize the major changes here. 
This work represents the first time a burning model which 
correctl y reproduces the nuclear statistical quasi-equilibrium 
(NSQE, iKhokhlovl [19891) to NSE transition time scales and 
length scales during incomplete silicon burning has ever been 
used in a multi-dimensional SN la calculation, as validated 
by comparison to steady-state detonation structures calcu- 
late d out to the pathologica l point with the ZND model (see 
e.g. iFickett & DavisI 119791) . Accurate reproduction of this 
low-density burning regime is essential because a signifi- 
cant portion_ofthe_^^Ni is produced in incomplete burn- 
ing dBravo et al.ll2010 ) so that the overall ^^Ni yield is de- 
termined by the details of how this low-density cutoff of ^^Ni 
production occurs. 

Obtaining a satisfactory reproduction of ZND detonation 
structures involved two main changes to the 3-stage burn- 
ing model. First, we found that the progress variable repre- 
senting the NSQE to NSE transition, 0^„, which also gives 
the mass fraction of Fe-group (NSE) material, did not match 
the time and space structure of this transition in steady state 
detonations calculated with a full nuc lear network. Th e ki- 
netics for this stage, first proposed by 'Khok hlovl (11991 ) and 
adopted in Calder et al. (2007), are given by the simple form 
d4>qn/dt = (1 - (j)qn) /tnse whcrc tnse is a calibrated timescale 
that depends mainly on temperature. We have found, how- 
ever that a far more appropriate match to steady-state det- 
onations at densities important for incomplete silicon burn- 
ing, p < 10^ g cm"-', is obtained with the alternative form 
d(j)q„/dt = (I - 4>gn)^ / TNiE- This necessitates a recalibration of 
Tnse, since it now plays a different role, but it is still sufficient 
for it to depend only on temperature. Although this change in 
derivative significantly improves the match between the pa- 
rameterized burning and the ZND structure at the densities 
of interest, it is still not exact at all densities. It is therefore 
unclear if this form is indicative of some underlying physi- 
cal process, and whether or not it is specific to detonations. 
There are several relaxation processes proceeding simultane- 
ously, so that it is non-trivial to quantify separate contribu- 
tions. This will be investigated in more detail in future work 
on post-processing nucleosynthesis. 

The second major burning model change was motivated by 
needing to match the thermodynamic, i.e. T, p, profiles at 
densities at which the portion of the detonation structure rep- 
resenting the NSQE to NSE transition is resolved on the grid. 
Our previous treatment released all of the nuclear binding en- 
ergy change to the NSE state by the end of the second stage, 
leaving the third stage NSQE to NSE transition energetically 
inert. This lack of an energy release on the NSE timescale 
leads to an incorrect progression of the density fall-off be- 
hind the shock front in the detonation. The abbreviated en- 
ergy release leads to an under-prediction of the temperature 
immediately behind the unresolved earlier burning stages in 
a propagating detonation. A very good match of thermody- 
namic profiles to the full-network steady state detonation was 
obtained by tying the completion of the energy release to 0^„, 
so that energy release occurs in three distinct stages. Note that 
the previous change involving the kinetics used for (j)q„ is also 
an important contributor to the realism of the thermodynamic 
profiles obtained. 

Finally, although detailed nucleosynthesis based on post- 
processing tracer particle histories will be published in fu- 
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ture work (Townsley et al., in prep.), we have performed 
two important tests for the burning model used in this work. 
The composition profiles obtained from post-processed his- 
tories for hydrodynamic simulations of steady-state detona- 
tions match the steady-state structure calculated via the ZND 
method with a large network with remarkable accuracy. Ad- 
ditionally, a preliminary version of the post-processing under 
development has been applied to the calculations presented 
here, and we have found good agreement between 0^„ and 
the the fraction of material in the form of Fe-group nuclides 
found from post-processing on both an overall basis and in 
ejection velocity bins. Overall, the burning model changes 
led to a modest but significant (around Q.2Mq) increase in the 
overall Fe-group yields for the same explosion. The yields 
found here are not q uite this much higher than similar cases 
from Townsley et al.l(2009) because the progenitors used here 
have a lower central carbon fraction. 

3.2. Improved Detonation Ignition Conditions 

In order to study the systematic effects associated with 
changing the DDT density, we need to minimize any sys- 
te matics in our m e thod o f starting a detonation. Previously, 
in iTownsley et al.l (l2009l) . we visualized the simulation data 
from the deflagration phase and plotted a density contour at 
the DDT density. When the flame reached k, 64 km away 
from the contour, we picked a computational cell half-way 
between the flame front and the DDT density contour to place 
a detonation ignition point with a radius of 8 km. The place- 
ment and size of the detonation point was chosen to be as 
close to the flame front as possible while still allowing the 
detonation point to develop into a self-sustained, stable det- 
onation front. If the detonation point is placed too close to 
the flame front, then the flame will interact with the detona- 
tion point before it develops into a self-sustained detonation. 
Comparisons of simulations from identical initial conditions 
with 8 and 12 km ignition radii finds the NSE yield differs 
by ^ 0.5% throughout the evolution. This results indicates 
that the total yield is insensitive to the choice of detonation 
ignition radius for radii less than the characteristic size of a 
rising plume (as can be seen in Figure|2]l. We adopt 12 km for 
the detonation ignition radius in our study as it produces more 
robust detonations at low density. 

To ensure that we do not introduce unintended systematic 
effects in this study, we improve our method of detonation ig- 
nition point placement over the previous "by hand" method 
by precisely defining the criteria for a DDT that is used in 
an algorithm. Parameters in th is algorithm are chosen to be 
consistent with iTownslev et aTl (12009.) . Once the flame front 
reaches the specified DDT density in a simulation, a detona- 
tion is ignited 32 km radially outward away from the flame 
front. Here the reaction-diffusion (RD) front is defined by the 
variable representing progress of the subsonic burning wave, 
0RD, with the leading edge defined as the region between the 
values 0.001 and 0.01 of this variable. During the deflagration 
phase of a simulation, (/ird is equivalent to the carbon-burning 
reaction progress variable, (fifg. At the leading edge of the RD 
front, very little carbon has burned and the local density is 
approximately equal to the unburned density. This provides 
a definition of the DDT density that is much more precise 
and accurate. We chose these criteria that ignite the detona- 
tion ahead of the RD front to avoid any issues with the det- 
onation ignition point overlapping with the artificially thick- 
ened flame. If we were to choose criteria that would initiate a 
detonation inside a thickened flame, the detonation structure 



would need to be altered in some way to be consistent with 
the artificial nature of that region. 

Our detonation ignition conditions also restrict detonation 
ignition points to be at least 200 km away from each other. 
This choice ensures that each rising plume starts 2^3 deto- 
nations, which is consistent with Townslev et al.l (120091) . In 
the case that multiple points within 200 km meet the deto- 
nation ignition conditions, the point furthest from the center 
of the star is preferred ensuring the ignition of detonations 
on plume "tops". In reality, the location of the spontaneous 
detonation points is not wel l known and i s the s ubject of ac- 
tive research. For instance, iRopke et al.l (120071) argue that a 
spontaneous detonation is triggered by the extreme turbulence 
found in the roiling fuel underneath the plume caps. Under the 
assumption that the DDT occurs when the flame reaches the 
low density for distributed burning, we place ignition points 
on the tops of the rising plumes. Future studies, however, 
will explore other physically motivated detonation methods 
in which the location of the detonation ignition is not neces- 
sarily specified relative to a plume, but rather determined by 
the local turbulence field, composition, dens ity, and tempera- 
ture (Ropke et al. 2007; Wooslev et al.ll2009h . Figure |2] shows 
the placement of first one detonation point, and then another 
more than 200 km away from the first, both at the specified 
DDT density. 
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Fig. 2. — Snapshots of realization 2 just after one (left) and then another 
(right) detonation was ignited at the specified transition density of 1.26 X 10 
g cm"^ (green contour). The reaction progress variable representing carbon 
burning is in color The detonating regions are the rapidly expanding regions 
of completely burned carbon ahead of the plume. 



Each detonation ignition point is defined by the profile of 
the reaction progress variable representing carbon burning, 
4) fa- Within the detonation ignition radius 0y„ = 1, and, be- 
cause we are igniting in fuel, (/>/„ = outside this region. The 
subsequent energy release from this change in 0/„ drives a 
shock strong enough to ignite the surrounding unburned ma- 
terial. This top-hat profile is a simple approach to igniting 
a detonation and has some drawbacks. At densities below 
« 10^ g cm"^, this detonation ignition method does not pro- 
duce a sufficiently strong shock to burn the fuel. The im- 
provements to our DDT method place the detonation point ra- 
dially outward from the specified DDT density; therefore, the 
density within the detonation point is actually somewhat less 
than that specified b y the DDT density. Th is differs from the 
method described in lTownslev et al.l (l2009h in that the detona- 
tion point was placed at a density that was somewhat higher 
than the specified DDT density. Igniting a detonation speci- 
fying a DDT density below 10^'^ g cm"^ is impossible using 
a top-hat profile because the detonation ignition point is actu- 
ally placed K.Kf g cm"^. In future works, we will explore the 
use of a gradient in 4>fa motivated bv lSeitenzahl et al1(l2009ah 
to describe the detonation ignition point profile which should 
result in the formation of stronger shocks that wifl burn the 



surrounding low-density fuel. 

4. PROPERTIES OF STATISTICAL SAMPLE AND METHOD 

To Study the systematic effect of transition density on the 
^^Ni yield synthesized during a simulated ex plosion, we uti- 
lize th e statistical framework developed in iTownslev et alj 
(|2009|). In ord e r to c ompare results between this study and 
iTownslev et alj (12009) . we use the same sample population 
using the same initial seed. The initial seed defines the starting 
point to a stream of random numbers used to characterize our 
sa mple population of the rmonuclear supernovae as described 
in ITownslev et al.l (|2009). In that study, we found that the 
sample dispersion in the estimated NSE yield does not asymp- 
tote until more than 10 realizations. Accordingly, our sample 
is made up of 30 realizations. 

During the early part of the deflagration phase (« 0. 1 s), the 
flame is most affected by convective motion in the core of the 
WD. Because the velocity field in the core and the number of 
ignition points are largely uncertain, we choose to character- 
ize the initial flame surface using spherical harmonics, each 
with a random coefficient picked sequentially from the initial 
seed. Each realization is defined with a unique perturbation 
on the initial spherically symmetric flame surface using 



a non-negUgible amount and we can more easily analyze the 
dependence of the yield on DDT density. 



'=/min 



(3) 



where ro is the radius of flame surface, gq is the amplitude of 
the perturbations. A; is a randomly-chosen coefficient corre- 
sponding to the spherical harmonic Y^, and /min and /^ax set 
the range of spherical modes used to perturb the flame sur- 
face. This method serves to initialize Rayleigh-Taylor unsta- 
ble plumes of random relative strengths, which we would ex- 
pect from various distributions of ignition points and varying 
strengths of the convective velocity field found in the real pop- 
ulation of progenitors. In this study, ro = 150 km, ao = 30 km, 
/mill = 12, and /^ax = 16. The choices of these parameters 
are motivated by the resolution of our study and the desire 
to obtain reasonable ^^Ni yields inferred fro m observations. 
These choices are discussed in further detail in lTownslev et al.l 
(l2009h . 

We choose to analyze the dependence on transition den- 
sity by choosing 5 different transition densities equidistant 
in log space at Cpooj = {7.1,7.2,7.3,7.4,7.5} in cgs units, 
where CpooT = logio Pddj- The range of DDT densities cho- 
sen for this study are motivated both by previous work and 
computational challenges. Below Cpooj = 7.1, a more re- 
alistic detonation structure is needed to successfully launch 
a detonation wave as desc ribed in the previous section and 
in ISeitenzahl et al.l (l2009ah . While recent studies have sug- 
gested a wider range of DDT densities are possible (iPan et al.l 
|T008; Schmidt et al. 2010), we suggest that trends resulting 
from varying the DDT density are captured with a maximum 

'CpDDT = 7.5. 

A simulation is performed for each of the 30 realizations 
in our sample at each DDT density for a total of 150 simu- 
lations. We choose not to explore DDT densities below this 
specified range for computational reasons. Because of the ap- 
proximate power-law density profile of the WD, DDT den- 
sities were chosen in log space because these densities cor- 
respond to relatively evenly spaced radial coordinates of the 
WD. The amount of ^^Ni synthesized in the explosion prin- 
cipally depends on the amount of expansion during the defla- 
gration phase. Therefore, spatially separated detonation ig- 
nition conditions will allow the amount of expansion to vary 
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Fig. 3. — The distributions of the ^*Ni yield for three different transition 
densities pv = 0.76 (solid), 1.26 (dashed), and 2.00 (dotted). Note that the 
variance increases with decreasing transition density. The distribution con'e- 
sponding to p^ = 0.76 is calculated by extrapolating the dependence of DDT 
density as described in JiS.ll such that the mean M(^*Ni) ~ 0.60 Mq consis- 
tent with observations. 

We find that extrapolating our results to a DDT density 
of CpDDT = 6.83 yields an estimated average ^^Ni yield of 
~ 0.60 Mq with a standard deviation of 0.21 Mq (from 
the fits described below and listed in Table IH l that is con- 
sistent with observations (JHowell et al.ll2009h . Because the 
actual DDT de nsity is unknown and the subject of on- 
going research (Pan et al. '2008'; 'AspdenetalJ l2008l I20T0I: 
IWooslev et al] [2009 ; Schmidt et a l. 20 13), we are free to 
choose this value of the DDT density as the fiducial DDT 
density, ^Pddt.o = 6.83. This choice is relevant for analy- 
sis and comparison to other works as discussed in ^ The 
distribution of ^''Ni material synthesized during the explo- 
sion is shown in Figure [3] for different transition densities at 
>CpDDT = 6.83, 7.10, and 7.30. 

The NSE mass, Mnse, is defined as J (j>g„pdV integrated 
over the star. We determine the NSE yield by running the 
simulation until Mnse is no longer increasing as a function of 
simulation time. This condition is defined as 



^^NSE ^QQ^ Mq 

dt s 



(4) 



Because we estimate the ^^Ni yield as a fraction of the NSE 
yield, we consider the ^^Ni yield to have plateaued when the 
NSE yield has plateaued. Considering additional ^''Ni from 
incomplete Si-burning in NSQE material and the efficient 
capture of excess neutrons onto Fe-group elements changes 
the final ^''Ni estimate by < 1%. Therefore, to good approxi- 
mation, the ^^Ni yield is a fraction of the NSE yield. Figure|4] 
shows the evolution of the NSE yield as a function of simu- 
lation time showing the DDT time and the NSE yield plateau 
time for each DDT density for realization 2. Discernible from 
this figure, there is a clear dependence of the NSE yield on 
transition density. 

5. RESULTS 

Within the DDT paradigm for thermonuclear supernovae, 
the conditions under w hich a transition occ urs are still not 
completely understood (JAspden et al.l 1200^ 2010). Gener- 
ally, though, a hypothesized transition occurs when the flame 
enters the distributed burning regime, which occurs when 
the flame speed equals the turbulent velocity at the scale 
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Fig. 4. — The NSE yields for realization 2 at each transition density used 
in this study. The closed circles indicate the time the flame first reaches the 
DDT density. The open squares show the time the NSE yield has plateaued 
as defined by Equation 0. 

of th e flame width assuming a K olmogorov turbulent cas- 
cade (iNiemever & WoosIevlll997l) . Recent work has placed 
more stringent requ i rements on the DDT dPan et alj 120081: 
IWooslev et all 120091: ISchmidt et ai] 1201 Oh and fundamental 
questions concerni ng deflagration in the hm it of disruptive 
turbulence remain (Polud nenko & Oranll2010 ). Regardless of 
the actual DDT mechanism, it is likely that the DDT condi- 
tions will depend on composition because both the width and 
burning rate of the flame depend on the abundances '^C and 
^^Ne. For the current study, we assume that a DDT will occur 
at a unique density given a particular composition, regardless 
of the microphysics involved. This assumption is reasonable 
given that the characteristics of the flame depend strongly on 
density. Using this assumption, we can delay the analysis of 
the particular microphysics that lead to a specific transition 
density and analyze the dependence of the amount of mate- 
rial synthesized to NSE during the explosion on the transition 
density. Therefore, we can think of each transition density 
as a proxy for changing the composition that determines the 
conditions for a DDT via t he appropriate microph ysics. 

As discussed in detail in lTownsley et alJ (l2009h . many sys- 
tematic effects exist that influence the outcome of a SN. In 
that study, the direct effect of ^^Ne was explored and found to 
have a negligible influence on the NSE yield. Therefore, we 
do not vary the initial ^^Ne mass fraction, but rather study the 
effect of varying the DDT density with the expectation that 
the NSE yield will be influenced indirectly by Ne and '^C 
abundances through the DDT density. We focus on the indi- 
rect effect of the ^^Ne abundance on t h e NSE yield, following 
up previous work in iTownslev et alj (l2009l) and neglect the 
effect of varying the carbon abundance. Additionally, we ne- 
glect effects due to the central ignition density, compositional 
and thermodynamic WD structure, and the total WD mass. 
The effects due to these variables will be studied in turn in fu- 
ture works with the goal of addressing interdependencies once 
individual effects are better understood. 

5.1. Dependence on Transition Density 

The evolution of the amount of material above a den- 
sity threshold (Mp-^p^^^^^ = J pdV, where we take pthie.s = 

2 X 10^ g cm~^) principally determines the dependence of the 
NSE yield on DDT density. This is due largely to the linear re- 
lationship between the NSE yield and Mp,>2(r = Tddt) shown 
in Figure |5] where fooT is defined as the time the flame first 
reaches the specified DDT density. This definition is consis- 
tent with our assumption that DDT conditions are met on the 
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Fig. 5. — Mass burned to NSE as compared to the mass above a density 
of 2 X lO' g cm"' at the first DDT time for each realization in our sample. 
The different shapes plotted (+, X, *, □, ■) correspond to different transi- 
tion densities (lO", lO'^, lO'-^ lO'", lO'^), respectively. The solid line 



shows the linear fit to yields produced with a DDT density of 10 g cm 



n-3 

most closely matching the DDT density used in Townslev et al. (2009). The 
dashed line shows a 1 : 1 coiTelation between the NSE yield and mass above 
2 X 10^ g cm"'. The lowest two DDT densities are less than the density 
threshold 2 X lO' g cm"' and thus show more scatter about the linear relation 
due to increased dependence on the plume morphology. 

tops of rising plumes. As mentioned above, other DDT loca- 
tions are possible, such as the highly turbulent region under- 
neath plume "caps", which would result in a higher DDT den- 
sity correlated to our present definition via the density struc- 
ture local to the rising plume. As shown in Figure H] the de- 
flagration phase burns only a relatively small fraction of the 
white dwarf and the majority of material is burned during the 
detonation. Once a detonation has started, the propagation 
speed of the burning wave is much greater than the rate of ex- 
pansion; therefore, the NSE yield is essentially independent 
of the evolution of Mp,>2 for t > foDT- The number of deto- 
nation points and their corresponding distribution in time and 
location contributes to the variance of the relation between the 
NSE yield and Mp,>2. The relation between the NSE yield 
and DDT density via Mpj>2 also depends on the acceleration 
of the RT-unstable plumes not being too great near foDT- For 
a constant plume rise-rate near Tddt, there is a linear relation- 
ship between the DDT density and foDT- This relationship 
allows the evolution of Mp,>2 to translate directly into a de- 
pendence of the NSE yield on DDT density. Figure |6] shows 
the evolution of Mp^y2 and pmin as a function of simulation 
time for 3 different realizations where pmin is the minimum 
density at which the flame is burning material. The evolu- 
tion of pmin shows the linear relationship between the log of 
the DDT density and fooT where the filled circles represent 
the foDT times for each of 5 transition densities and the filled 
squares show the fiducial DDT density. For the times corre- 
sponding to the DDT conditions, the evolution of Mp,>2 is 
in a region where the rate of expansion of the star becomes 
significant and Mp,>2 begins to drop off relatively quickly. 

We find, for each individual realization, the NSE yield de- 
pends quadratically on transition density. Table [3] lists the 
NSE yields for each realization at each of 5 transition den- 
sities equidistant in log-space. Table |4] shows the coefficients 
used to fit the quadratic dependence of NSE yield on transi- 
tion density for 



M^se(Mq) = a (Cpddt) + bCpooT + c. 



(5) 



Figure |2] shows fits to the NSE yield for each realization and 
the average yields at each transition density. The NSE yields 
for individual realizations are plotted as blue crosses and the 
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Fig. 6. — Comparison of the evolution of the white dwarf defined by the mass above the density threshold (Mp,>2) to the plume evolution defined by the 
minimum flame density (pmin)- The left panel presents the evolution ofMp^^2 (solid) and pini„ (dotted) as a function of simulation time for realizations 2 (blue), 
10 (red), and 18 (green) showing the expansion rate of the WD and the plume rise-time with DDT times (circles) associated with each DDT density emphasized. 
Mpj>2 is defined as the total mass with density greater than 2 X 10^ g cm"' and pn,in is defined as the minimum density at which the flame is burning material. 
The translation into Mp^^i as a function of p„,in is shown in the right panel for comparison to A/^se as a function of poDT shown in Figure|2] We also highlight 
the fiducial DDT density as black squares. These plots emphasize the importance of the rate of expansion of the WD and the plume rise-time as a function of 
density in determining the relation between the NSE yield and DDT density. Realization 2 produces NSE yields that are close to the average while realizations 
10 and 18 produce the lowest and highest yields, respectively. 



individual curves showing the dependence of NSE yield on 
DDT density are shown in grey. The average NSE yields at 
each DDT density are shown in red. At a transition density 
of 10^^ g cm"^, notice that the dependence on DDT density 
has flattened out for almost all realizations and the variance 
of the yields among all realizations is small. This behavior 
is due to the fact that the progenitor white dwarf star has a 
finite mass of 1.37 Mq. At the highest transition density, the 
expansion rate of the star is relatively small, and the evolution 
of Mp,>2 is relatively slowly declining as compared to lower 
DDT densities (see Figure|5ll. 



TABLE 3 

NSE YIELDS IN Mq FOR EACH REALIZATION AT EACH TRANSITION 

DENSITY, £pDDT- 
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Fig. 7. — Quadratic fits (grey lines) to the NSE yields for each realiza- 
tion (blue crosses) and the average NSE yield at each transition density (red 
circles). The qualitative similarity between this figure and the right panel of 
Figure|6]is due to the near 1 : 1 correlation between Mnse and Mp^^2(f = 'ddt)- 



The curvature of the NSE yield dependence on DDT den- 
sity, a, is well correlated with the NSE yield at a given DDT 
density. Figure |8] shows this correlation for CpoDT.o = 6.83; 
however, a correlation exists for all DDT densities as may be 
discernible from Figure |7] noting that most black lines rep- 
resenting individual realizations do not cross. The lower the 
NSE yield for a given realization, the stronger the dependence 
on transition density. This result is likely due to realizations 
with lower yields having multiple competing plumes that re- 
lease more energy allowing the star to expand more rapidly 
leading to a stronger dependence on the transition density. 
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7.1 


7.2 


7.3 


7.4 


7.5 


1 


0.989 


1.082 


1.152 


1.207 


1.236 


2 


0.970 


1.056 


1.131 


1.184 


1.215 


3 


0.948 


1.043 


1.119 


1.174 


1.209 


4 


0.951 


1.062 


1.144 


1.197 


1.226 


5 


1.027 


1.101 


1.168 


1.214 


1.242 


6 


1.074 


1.138 


1.185 


1.219 


1.248 


7 


1.076 


1.132 


1.179 


1.215 


1.243 


8 


1.110 


1.167 


1.210 


1.239 


1.260 


9 


0.917 


1.022 


1.105 


1.163 


1.199 


10 


0.791 


0.929 


1.038 


1.113 


1.166 


11 


1.092 


1.147 


1.189 


1.229 


1.255 


12 


1.240 


1.258 


1.273 


1.289 


1.301 


13 


0.882 


0.987 


1.083 


1.145 


1.183 


14 


1.070 


1.125 


1.185 


1.223 


1.250 


15 


1.119 


1.176 


1.218 


1.254 


1.274 


16 


0.883 


0.994 


1.081 


1.142 


1.183 


17 


1.026 


1.115 


1.171 


1.213 


1.241 


18 


1.181 


1.223 


1.249 


1.270 


1.290 


19 


0.967 


1.057 


1.132 


1.185 


1.220 


20 


1.005 


1.089 


1.138 


1.180 


1.218 


21 


1.222 


1.248 


1.268 


1.284 


1.297 


22 


1.045 


1.104 


1.170 


1.216 


1.244 


23 


1.001 


1.085 


1.150 


1.187 


1.223 


24 


0.900 


1.012 


1.091 


1.130 


1.176 


25 


0.989 


1.057 


1.139 


1.193 


1.226 


26 


1.092 


1.150 


1.192 


1.228 


1.254 


27 


1.192 


1.222 


1.251 


1.274 


1.290 


28 


1.145 


1.183 


1.220 


1.253 


1.273 


29 


1.063 


1.131 


1.175 


1.210 


1.235 


30 


1.179 


1.225 


1.249 


1.273 


1.288 


M 


1.038 


1.111 


1.169 


1.210 


1.239 


(T 


0.110 


0.082 


0.059 


0.046 


0.036 


'^M 


0.020 


0.015 


0.011 


0.008 


0.007 
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TABLE 4 

Coefficients used to fit the quadratic dependence of 

transition density on nse yield using equation (o. the nse 

MASS IN UNITS OF Mq IS EVALUTAED AT £pDDT,0 USING THE 
COEFFICIENTS. 



b = Sc + j, 



(7) 



# 


a(MQ) 


b(MQ) 


c(Mq) 


MnseC^Pddt.o) 


1 


-1.024 


15.56 


-57.9 


0.637 


2 


-0.933 


14.24 


-53.1 


0.633 


3 


-1.010 


15.40 


-57.5 


0.589 


4 


-1.388 


20.95 


-77.8 


0.515 


5 


-0.815 


12.43 


-46.2 


0.732 


6 


-0.595 


9.12 


-33.7 


0.852 


7 


-0.480 


7.42 


-27.4 


0.876 


8 


-0.619 


9.41 


-34.5 


0.898 


9 


-1.164 


17.70 


-66 A 


0.515 


10 


-1.454 


22.17 


-83.3 


0.279 


11 


-0.442 


6.87 


-25.4 


0.902 


12 


-0.083 


1.36 


-4.3 


1.183 


13 


-1.212 


18.45 


-69.0 


0.456 


14 


-0.557 


8.58 


-31.8 


0.844 


15 


-0.575 


8.78 


-32.2 


0.910 


16 


-1.173 


17.88 


-66.9 


0.470 


17 


-0.963 


14.58 


-54.0 


0.712 


18 


-0.348 


5.34 


-19.2 


1.048 


19 


-0.926 


14.15 


-52.8 


0.628 


20 


-0.708 


10.86 


-40.4 


0.740 


21 


-0.210 


3.25 


-11.3 


1.134 


22 


-0.583 


9.03 


-33.6 


0.798 


23 


-1.093 


16.47 


-60.9 


0.696 


24 


-1.219 


18.46 


-68.7 


0.502 


25 


-0.782 


12.04 


-45.1 


0.671 


26 


-0.490 


7.56 


-27.9 


0.896 


27 


-0.248 


3.86 


-13.7 


1.080 


28 


-0.276 


4.36 


-15.9 


1.006 


29 


-0.685 


10.42 


-38.4 


0.827 


30 


-0.429 


6.54 


-23.6 


1.032 
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Fig. 8. — The correlation between the NSE yield at the fiducial DDT den- 
sity CpDDio = 6.83 and the fitting parameter a from Equation (5). A lower 
yield indicates a stronger curvature with DDT density. 

After fitting the dependence on DDT density for each real- 
ization, we find an interesting correlation between the fitting 
parameters from Equation Q. Figure |9] shows the correlation 
between the fitting parameters a, b, and c with c on the hor- 
izontal axis. The tight correlation between these parameters 
indicates that a single parameter describes the dependence on 
DDT density for a given realization. The fitting parameters a 
and b can be expressed as a function of c: 



where a = 1.754 ± 0.003 x lO'^, ^ = -0.005 ± 0.002, S = 
-0.2646 ± 0.0005, and 7 = 0.21 ± 0.02. These parameters 
were calculated using a Least Squares method and the associ- 
ated errors were calculated from the corresponding covariance 
matrix. 
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Fig. 9. — The correlation between the NSE yield fitting parameter, c. with 
the fitting parameters a (left axis) and b (right axis) from Table|4]for each re- 
alization. The fitting parameters are defined by Equation J5) and the relations 
between the fitting parameters are given by Equations J6j and (7). 



Each realization has a different random initial perturbation 
of the central ignition condition that sets the plume morphol- 
ogy. Dominant single plumes tend to allow the star to expand 
less prior to reaching the conditions for a detonation, while 
multiple competing plumes tend to release more energy dur- 
ing the deflagration phase, allowing more expansion prior to 
reaching the conditions for a DDT. Shown in Figure [TO] are 
the Rayleigh-Taylor unstable plumes for the realizations with 
the highest and lowest yields a few tenths of seconds into the 
deflagration. Unfortunately, no strong correlation exists be- 
tween the properties of the initial flame surface, such as mass 
enclosed, surface to volume ratio, or amount of power in the 
perturbation, for a given realization and the single fitting pa- 
rameter c. This lack of a correlation implies there is no way 
to tell whether a particular initial condition will seed a single 
dominant plume or multiple competing plumes. 

Physically, we might expect two competing effects that in- 
fluence the NSE yield: plume morphology (rise time) and the 
rate of expansion. In our 2D simulations these two effects ap- 
pear to be correlated as seen in Figure |6] Our results indicate 
that plumes near the symmetry axis tend to rise faster than 
plumes rising near the equator and we attribute this result to 
the fact that our simulations are two-dimensional. A plume 
that develops near the symmetry axis naturally becomes dom- 
inant in 2D and determines the DDT time and thus the NSE 
yield. In addition, a plume near the symmetry axis represents 
less volume than a plume of similar size near the equator This 
indicates that energy is being deposited into a smaller volume 
allowing for a lesser rate of expansion. This explains why the 
overall rate of expansion of the star is correlated to the rise 
time of the first plume to reach DDT conditions. 

In order to evaluate the dependence on DDT density with- 
out this unphysical correlation, we need to perform 3D sim- 
ulations. A suite of 3D simulations will likely create a two- 
parameter family (plume rise and expansion rates) of solu- 
tions which describe the yield as a function DDT density. Ad- 
ditionally, we expect 3D plumes to behave similarly to 2D 
on-axis plumes implying a faster plume rise time. 
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Fig. 10. — Snapshots of realization 18 (left) and realization 10 (right) at 
a simulation time of 0.4 seconds. Reahzation 18 produced the highest yield 
and shows the development of a single dominant plume while realization 10 
had the lowest yield and shows all plumes developing at about the same rate. 
Shown in color are fuel and burning products: unburned C, O, Ne (yellow) 
and Fe-group (NSE, black). Density in g cm"' is indicated by contours (blue) 
logarithmically spaced at integer powers of 10 as well as the DDT density of 
1.26 X 10''gcm"3 (red). 



TABLE 5 
Statistical properties of the fitting parameter c and the NSE 

YIELD at PdDT = 6.76 X 10* G CM"'. 



c(Mq) 



Mean 



-42 
0.77 



Std. Dev. 



21 
0.22 



Std. Dev. of Mean 



4 
0.04 



As a result of the tight correlation for a and b and the de- 
termination of Equations (1617b . the average relation is char- 
acterized by just the average c — or equivalently, the average 
NSE yield at a specified pddt- We evaluate the average and 
standard deviation of the fitting parameter c as well as the sta- 
tistical properties of the NSE yield at the fiducial transition 
density of 10^^^ g cm"^ shown in Table|5] 

5.2. NSE Yield Dependence on ^^Ne 

Recall that for the purposes of this study, we are neglecting 
the effects of varying the core carbon abundance, the central 
ignition density, the WD structure, etc. except for the indirect 
effect of X22, the neutron excess produced during simmering, 
and the DDT density. In order to derive the NSE yield depen- 
dence on X22, we need to expand the derivative of Mnse with 
respect to X22 which involves only a couple of terms given our 
assumptions 

^Mnse 



9Mnse dCpxim: dM^sE 



(8) 



dX22 dCpuDT dX22 9X22 

iTownslev et all (l2009h showed that ^^e = o. In that study, 
they employed an estimate of the mass burned to NSE by mea- 
suring the amount of mass above a density of 2 x 10^ g cm""*. 
That correlation is confirmed by our current study as shown 
in Figure |5] where we plot the NSE yield of all realizations at 
each transition density. The correlation between NSE yield 
and mass above 2 x 10^ g cm"-' at foDX for a transition den- 
sity of 1.26 X 10'' g cm"-', closest to the transition density 
used for that study, is highlighted. Given our assumptions 
and the result that th e NSE yield does not depend directly on 
X^2 dTownslevet all 12009 ). the NSE yield only depends di- 
rectly on the DDT density. The affect of Xjj on the yield 



TABLE 6 

COEpnCIENTS FOR LOG-LOG FITS TO EQUATION (fTTt . THE SUBSCRIPT 

DENOTES WHETHER FITTING FOR THE LAMINAR FLAME SPEED (s) OR 

FLAME WIDTH (5). 





X22 = 0.02 


X22 = 0.06 


as 


0.7942 


0.7745 


bs 


-12.735 


-12.121 


as 


-1.3550 


-1.3507 


bs 


19.17 


18.88 



enters through the DDT density. Therefore, we will construct 
the functional dependence of Cpddt on Xjj. 

The DDT density depends on X22 via the microphysics in- 
volved in determining the conditions under which the flame 
transitions from a deflagration to a detonation, which, as noted 
above, is incompletely understood. We assume the transition 
occurs at a particular density at which the flame enters the dis- 
tributed burning regime, but more stringent conditions for the 
DDT include dependences on the turbulent cascade and the 
growth o f a critical mass of fuel with sufficiently strong tur- 
bulen ce dPan et alj|2008l: IWooslev et ai]|2009l: ISchmidt et all 
I2OIOI) . The dependence on composition of these models may 
be explored and applied to the trends with DDT density pre- 
sented in this study. Under our assumptions, the flame enters 
the distributed burning regime when the laminar flame width 
becomes of order the Gibson length (/g) where 



\u' 



(9) 



Here, L is the length scale on which the strength of the tur- 
bulent velocity (u') is evaluated, and si is the laminar flame 
speed (see e.g. iPetersI 12000.) . By using th e laminar flame 
speeds and widths which depend on X22 from lChamulak et al] 
( l2007h . we choose the conditions for a detonation are met 
at X22 = 0.02 at a particular transition density {CpBui) when 
the Gibson length is equal to the laminar flame width. The 
strength of the turbulent velocity field is calculated assuming 
a Kolmogorov turbulence cascade: 



u' = si{p,X22 = om) 



5,{p,X22 = Qm) 



1/3 



(10) 



Using m', we solve for the change in transition density, 
A£/3ddt, by changing ^^Ne to X22 = 0.06 and again setting 
the Gibson length equal to the laminar flame width where L 
cancels out of the equation. This evaluation yields a change 
in transition density over a change in X22, or J^°'' . 

We use log-log fits to the flame speed and flame width 
as a function of den sity using the table generated by 
IChamulak et al.l (l2007h . We use only densities above 10** g 
cm"3 for Xn = 0.3 and above 5 x 10** for Xn = 0.5. Within 
this parameter space, a power-law dependence of the flame 
speeds and widths on the log of the density is well-defined 
and the algorithm used to solve the flame characteristics is 
more stable at higher densities. Our resulting expressions for 
the density dependence on flame speed, s, and flame width, 
5i, at the carbon mass fraction used in this study {X\2 = 0.4) 
for X22 = 0.02. 0.06 are given by 



\vis = a\vip+b 
with coefficients given in Tabled 



(11) 
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Our derived expression for the derivative of transition den- 
sity as a function of X22 is given by 

dCpDDT bsfi-bs.i + 'i {bs.i-b,,6)+\\\pr,Y)T:{'ias.2-as^2) 



dX22 



In 10 (3fl.s.6 - flA-,6) AX22 



^p- 



'DDT 



AX22 
■■uCpDDT + V, 



(12) 



where u = 0.43 15 and v = -6.301 . We solve this first-order dif- 
ferential equation and express the DDT density as a function 

0fX^2. 

CpuDT(X^2)=l{e"^''''~'''''''^- 

+CpDDjfle"^''''-'''"'\ (13) 

where X22 is the parameterized ^^Ne mass fraction chosen to 
be associated with the fiducial transition density, Cpddj.o- For 
all the transition densities considered in this study, the inter- 
face between the core composition and the composition of the 
outer layer is at a lower density. Therefore, the relevant -^Ne 
content to consider is that of the core. Plugging in Equation 
( fT3] l into Equation (|5j, we obtain the functional dependence 
of Mnse onX22, such that 

Mnse = Mnse {Cpddt (^22) ) ■ (14) 

Equation (fT4l l is evaluated and plotted in Figure [TT] for the 
fiducial transition density at the ^^Ne mass fraction used in 
this study, 

CpBDT(X^2.0 = 0-03) = CpDDT^O = 6.83. (15) 

We propagate the standard deviation of the mean evaluated at 
CpDDj.o for a range of X22. This is calculated by considering 
the relation between the standard deviation of the mean of the 
NSE mass and the standard deviation of the mean of the fitting 
parameter, c, given by 

o-NSE = (a [CpDDT {X^2)f + SCpDm (^22) + 1) c^o (16) 

where unse is the standard deviation of the mean of the NSE 
mass and CTc is the standard deviation of the mean in the fitting 
parameter c. We evaluate CTc by inverting Equation ( fTSI l and 
solving for X22 = ^22 q. Plugging in this solution, the standard 
deviation of the mean of the NSE mass as a function of X22 
becomes 



cnse (^^22) - 



a [Cpum (X^2) ] " + ^^PDDT (^22) + 1 

a [>CpDDT,o] + SCpDDT.O + 1 



O'NSE,0- 



(17) 



5.3. 



^^Ni Yield Dependence on Metallicity 

Now that we have constructed the functional dependence of 
NSE yield on ^-Ne we need to consider the dependence of the 
amount of ^''Ni synthesized in the explosion on metallicity 
through the ^^Ne content. A fractional amount of the NSE 
material is radioactive ^''Ni, which powers the supernova light 
curve. This fraction depends on Yg (and ^^Ne) as described in 
lTimmesetan(l2003h : 



M56 = /56MnsE- 



(18) 



In order to determine the dependence of M^f, on metal- 
licity through the DDT density, we must construct /se and 
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22 

Fig. 1 1 . — The analytic solution to Mnse as a function of parameterized 
X!,2 and using standard error propagation to obtain the uncertainty based on 
the standard deviation of the mean (dashed lines). The vertical dot-dashed 
line indicates the parameter space in which this study was performed. These 
results were evaluated for the fiducial transition density of 6.76 X 10* g cm"^ 
corresponding to Zjj = 0.03. 

explore any dependencies on DDT density and metallicity. 
From our simulations, we calculate Mnse directly and esti- 
mate Mse from Y^ choosing the non-^'^Ni NSE material to be 
50/50 ^"^Fe and ^^Ni by mass. Recall that X22 from Equation 
([T]i contains a term, AX22, which is a parameterization of the 
change in Y^ due to neutronization during the carbon simmer- 
ing phase. However, not all of the neutron excess evaluated 
at the Mnse plateau time comes from ^^Ne or the carbon sim- 
mering products. Some change in Y^ is due to weak reactions 
occuri ng during the explosion that are included i n our burning 
model (ICalder et al.ll2007HTownsley et al ."2007'). 

First we consider the amount of non-^^Ni NSE material de- 
termined by X22 by equating the initial Y,, to the Y^, of material 
in NSE. Using baryon and lepton conservation for NSE ma- 
terial, we describe the electron fraction by contributions from 

54pg^ 58j^j ^ ^jjjj 56j^j^ 



F.= 



26 /I 



54 V 2 



r/ni 



28 /I 



58 V 2 



r/n( 



28 

+ 77(1 -/non-56), 

j6 



(19) 
where /non-56 is the mass fraction of non-^^Ni NSE material. 
For the following evaluation, we approximate the composition 
to be that of the core because most of the NSE material is 
within the core. We write the initial K as 
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Yp.i — ——Xnn + — ( 1 At 



^22 



22"- ■ 2 ^- --22) ■ (2^^ 

We can then solve for /non-56 by equating (fT9] l and (|20] | and 
add a term X„ to represent the additional contribution due to 
neutronization from weak reactions occuring during the ex- 
plosion 

/non-56 = T7j^^22 + ^« ■ (21) 



Then the Ni fraction of material in NSE is 



/56 = 1- 



783 
308 



Xoo Xf, 



^22 



(22) 



We note that the rate of weak reactions occuring during the 
explosion may depend on the initial composition as well. Ac- 
cordingly, X„ may have a dependence on X22. We also expect 
X„ to vary as a function of the transition density because a 
higher transition density will have less time for weak reac- 
tions to occur 

We construct a statistical sample of /ss using the ratio of 
the ^''Ni and NSE yields produced in the simulations. We 
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calculate the dependence of X„ on transition density using a 
least-squares method and the result is shown in Figure [121 
Evaluating the partial derivative of X„ at the fiducial transition 
density, we find a shallow dependence: 

= -0.096. (23) 

^Pddt=-Cpddt,o 



0.14 



dCp 



DDT 




0.08 



6.9 



7.1 7.2 
logloPoDT 



Fig. 12. — The fraction of non-^*Ni N SE m aterial due to weak processes 
(X„) is calculated from inverting Equation J22) and using the /sg = Msj/Mnse 
averaged over all realizations at each DDT density. The standard deviation 
of the sample (outer) and mean (inner) are shown as error bars at each DDT 
density. A quadratic best fit Une is calculated using the standard deviation 
of the sample. The fiducial DDT density is shown as a black circle where 
X„ = 0.133. 

We want to know the dependence of M56 on X22. Since 
/56 depends on the neutronization from weak reactions dur- 
ing the explosion, we must consider any dependencies of this 
neutronization on DDT density or X22. We show that the de- 
pendence of /56 on ^22 via the effect of X22 on the in-situ neu- 
tronization is much weaker than the direct dependence due 
to lepton number conservation. Similar to Equation ( fTSl l. we 
write 

dXii 9Mnse dX2i 8X22 ' 
Using Equation ( fTST l and expanding the partial derivative of 
M56 on X22, we obtain 



dMsf, c/Mnse ,^ dfsf, 

= /56 —TTr ^ -MnSE - 



(25) 



dX22 dX22 dX22 

Taking the derivative of Equation ( |22] | and expanding on X22, 
we obtain 

dfse _ dfs6 dX„ 

dX22 ~ aX22 dX22 ■ 
Now we wish to evaluate whether ^^ is an important contri- 
bution to the overall evaluation of the ^^Ni mass. Expanding 
this term yields 

dX„ 



dX„ dCpDDT dXn 



(27) 



^^22 dCpuDT dX22 9.X22 

Referring to iTownslev et al.l (l2009h . we can estimate -S^ by 
calculating the average ratio of M56 to Mnse at X22 = 0,0.02 
for the first 5 realizations whose detonation phases were sim- 
ulated. The result is ^^ ^ -0.2. The first term in Equa- 
tion ( |27] ) can be evaluated from multiplying Equation (1231 ) and 
Equation ( fT2b for CpuoT in the range 7.0-7.5. For the fiducial 
transition density (Cpddt.o)^ we find 
dX„ 



dX 



= (-0.126)(-3.35)-0.2-0.2. 



(28) 



Comparing to gj^ ~ -2.5, we find that the magnitude of -^ 
is much smaller and is unimportant for our evaluation of the 
dependence of the ^^Ni yield on X22. Therefore, we can ignore 



this term in the expansion of ^-M^ in Equation 



and let 



^f- ^ 0. The full derivative of /se with respect to X22 can 
now be written as a partial derivative such that 

^^=/,,f^^!^+M.sE|^. (29) 

dX22 dX22 0X22 

While we approximate Xn as constant, we choose to eval- 
uate it at iZpDDT.o using the best fit curve from Figure [12] ob- 
taining X„ = 0.133. Now w e can relate the meta llicity to ^^Ne 
since X22 traces metallicity dTimmes et alj2003l) . Substituting 
X22 = O.O14(Z/Z0) in Equation ([TJ, we obtain 



X^2 = 0.014 7— + AX22(AJ;), 



(30) 



where AX22 = 0.01 for this study and Xjj is our parametriza- 
tion of the actual ^^Ne mass fraction, X22. The ^^Ni yield as a 
function of X22 is calculated by multiplying Equation ( i22b by 
Equation ([T4l) and using X„ = 0.133. The ^^Ni yield is plotted 
in Figure [13] using Equation ( l30l l to relate X22 to Z/Zq. 




Fig. 13. — The solution of Mnse (red) and Mjg (green) computed as a 
functio n of metall icity as compared t o the '*Ni relations from Timmes et al] 
(20031) (blue) and IBravo et alj (20101) (magenta) normalized to the average 
'*Ni yield from our simulations. The dashed lines show the propagated stan- 
dard deviation of the mean. The vertical dot-dashed line indicates the param- 
eter space in which this study was performed. These results were evaluated 
with a fiducial transition density of 6.76 X 10^ g cm"^ at Z = 1 .4Zq . 



6. CONCLUSIONS AND FUTURE WORK 

We have analyzed the influence of the DDT density on 
the total ''^Ni synthesized during a thermonuclear super- 
novae. We determined that the dependence on DDT density 
is quadratic in nature, but this dependence can be described 
by a single parameter We estimated the dependence of the 
supernova brightness (^''Ni yield) on metallicity by assuming 
a DDT occurs when the flame enters the distributed burning 
regime and extrapolating the laminar flame speeds and widths 
down to low densities. We find 



dMib 



diZ/Zo) 



= -0.067 ±0.004Mp 



©J 



(31) 



Z=Zn, 



22 



which is slightly steeper than lTimmes et al.l (120031) as seen in 
Figure [T3] The uncertainty was calculated using the standard 
deviation of the mean of the fitting parameter c — or equiva- 
lently, the standard deviation of the mean of the NSE yield at 
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a particular DDT density. The uncertainty in the assumptions 
about the normalization of the transition density as a function 
of ^22 was not considered for the purpose of evaluating the 
uncertainty in the derivative. We find the effect of metallicity 
on DDT density influences the production of NSE material; 
however, the ratio of the ^^Ni yield to the overall NSE yield 
does not change as significantly and remains similar to the 
relation estimated from approximate lepton number conser- 
vation. We also find that the scatter in supernova brightness 
increases with decreasing transitio n dens ity. 

The very recent work of Bravo et al.| (|2010|) on metallic- 
ity as a source of dispersion in the luminosity- width relation- 
ship of bolometric light curves addresses many of the same 
issues as our stu dy and warrants discussion. In particular, 
iBravo et al] (1201 0) derive a metallicity dependence on the 
DDT density similarly to our study. However, the DDT den- 
sity in their ID simulations is very different from the DDT 
density in ours. The principal differenc e between the work 
described in this manuscript and that of iBravo et al.l (1201 Ol) 
follows from the use of multi-dimensional simulations. While 
3D simulations are required to correctly capture the effects of 
fluid dynamics and the Rayleigh-Taylor instability, 2D sim- 
ulations incorporate these effects and relax the assumption 
of symmetry. Breaking the symmetry and assuming that a 
DDT initiates as a rising plume approaches pddt produces 
an expansion history very different from what would be ob- 
served in ID simulations. In fact, using ID simulations im- 
plies that pddt plays an unphysical role. The real physical 
description of the DDT in a type la supernova will depend on 
flame-turbulence interactions that will themselves depend on 
multi-dimensional effects in the flow. ID models may be able 
to capture these effects, but such models must be motivated 
by more physical multi-dimensional studies. 

In addition, our statistical framework with realizations from 
randomized initial conditions allows calculation of the scat- 
ter inherent in the multi-dimensional models. The standard 
deviations calculated for our models averaged over a set of 
realizations demonstrate that there can be considerable varia- 
tion following from the randomized initial conditions. These 
variations follow from the different amounts of expansion oc- 
curring during the deflagration phase that follow from the dif- 
ferent rising plume morphologies. In ID models, even with 
the progenitor metallicity determining pddt, there is a one-to- 
one correspondence between expansion and pddt for a given 
proge nitor. 

Bra vo et all ( 1201 Oh explored two scenarios, a linear and a 
non-linear dependence of the ^^Ni yield on Z, and report ex- 
cellent agreement bet ween the non-linear s cenario and ob- 
servations reported by [Gallagher et al.l ( l2008h with the caveat 
that they used ID models to arrive at their conclusion. With- 
out compa ring results involv ing the metallicity dependence 
on pddt. Bravo et alj (12010) find a stronger dependence of 
the ^^Ni synthesized from incomplete Si-burning on Z than 
^^Ni coming from NSE material. We find a shallower depen- 
dence of the ^*'Ni yield on Z than ' Bravo et al.1 (1201 Oh in part 
due to our assumption that ^*'Ni is synthesized as a fraction 
of NSE material; however, comparing their stratified models 
to the Timmes et al. (2003 ) relation results in only a ^^ 20% 
difference in the dependence of ^^Ni yield on Z. This effect 
is compounded in their "non-linear" scenario, but again, we 
must emphasize that pddt plays a completely unphysical role 
in their simulations and the degree to which this effect is ac- 
tually enhanced is unknown. 

The light curve width calculations and subsequent popula- 



tion synthesis performed by iBravo et al] (1201 Oh after finding 
a dependence of the ^^Ni yield on metallicity is useful in es- 
timating a metallicity dependence in the bubble residual. For 
this purpose, w e com pare the results of the non-linear scenario 
of Br avo et aP ((2010) to our study along with the expec ted de- 
pende nce due to le pton number conse rvation (Timmes et alJ 
I2003h in Figure [H iBravo et all (1201 Oh find a - 30% steeper 
dependence of the ^^Ni yield on Z/Zq evaluated at Z = IZ© 
than our result in Equation dJTT l. In addition, they find that 
for high Z the ^''Ni yield tends to flatten out becoming more 
similar to our results. For a particular subrange of metal- 
licites, our results are very similar indicating t hat by per- 
forming the same analysis as iBravo et aP (1201 Oh our results 
should also agree with the metallicity dependence found by 
'Gallagher etal] JlOOSb. 

This similarity depends o n th e choice of rnean ^ ''Ni yield 
as both lTimmes et al.1 (l2003h and lBravo etall (1201 Oh relations 
are proportional to this quantity. This implies that the steep- 
ness of the relations are affected by the choice of mean ^*Ni 
yield. Our results are not sensitive to this choice. Choos- 
ing a higher fiducial DDT density of 10^ g cm"^ that results 
in a higher mean ^^Ni yield of 0.77 Mq increases Equation 
( [3TT i by only 0.005 Mq. We see less dependence on the mean 
^^Ni yield because, although we have a similar dependence 
on ^22 due to lepton number conservation, the dependence 
of the yield on Cpom is stronger at lower fiducial CpoDT.o, 
as shown in previous sections. This effect is not c aptured in a 
simple proportionality relation like that quoted by Bravo et al.l 
(1201 Oh . even if this effect is present in their calculations. 

It has been discussed whether the source of scatter in peak 
brightn ess itself can be attributed to metallic ity dTimmes et al.l 
2003; Howell et al.ll2009l; iBravo et alj|2010h . We submit that 
multi-dimensional effects following from fluid instabilities 
during the deflagration phase leading to varying amounts of 
expansion provide scatter consistent with observations. The 
influence of metallicity on the DDT density affects the du- 
ration of the deflagration which will secondarily influence the 
magnitude of the scatter Our results show that the primary pa- 
rameter is the degree of (pre-)expansion before DDT — which 
determines the amount of mass at high density. This is, in 
turn, controlled by both the expansion rate and the plume rise 
time. These are also expected to be the basic ingredients in 
reality. 

While our study is consistent with the expected theoret - 
ical brightness trend with metallicity dTimmes et alJ 120031) . 
observation s to date have not been able to confirm thi s 
prediction ( Gallagher et 10120051 120081: iHowellet alJ l2009'). 
The progenitor age is difficult to decouple from metallic- 
ity gi ven the mass-metal licity relationship within galaxy 
types (" Gallagher et al.l2008 ). Additionally, the dependence of 
mean brightness of SNe la on the age of the parent stellar pop- 
ulation app ears to be much stronge r than any depend ence on 
metaflicity dGallagher et al ...2008; Howell et al.ll2009h . Seem- 
ingly, the only way to observe a dependence on metallicity 
is to constrain the mean stellar age by selecting gala xies of 
the sa me type. This approach was used by Gallaghe r et alj 
(2008), but the difficulties in accurately measuring metallici- 
ties for the parent stellar population have only constrained the 
magnitude of the effect on mean brightness, but thus far have 
not proved that a metallicity effect exists. 

In order to determine the dependence of supernova bright- 
ness (^^Ni yield) on metallicity, a better understanding of 
how the transition density is affected by changes in the ^^Ne 
content is needed. Currently, we have extrapolated data 
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from I Chamulak et alj (l2007h down to the range of expected 
transition densities using flame dat a from densities above 
10** g cm"^ The trend observed by IChamulak et alJ (l2007h 
that increasing X22 increases the laminar flame speed is valid 
for densities above 10^ g cm"^; however, it is unclear if this 
trend will hold for densities below 10^ g cm"^. Direct nu- 
merical simulations of the flame are needed for these lower 
densities to determine whether this trend still holds. How- 
ever, in this parameter space, the laminar flame is extremely 
slow requiring a low-Mach-number treatment to model the 
flow. Regardless of the choice of model for the mechanism 
that produces a spontaneous detonation, a necessary condition 
is thought to be distributed burning. Properties of the lam- 
inar flame are necessary to estimate when burning becomes 
distributed. Future studies exploring other DDT mechanisms 
must also determine their compositional dependencies. So far, 
we have onl y explored the effect of varying the ^^Ne con- 
tent directly dTownslev et al.ll2009l) and indirectly through the 
DDT density. Other effects exist that we have yet to study 
such as the carbon composition and flame ignition density (set 
by the average progenitor age). These properties may be influ- 
enced by metallicity such that the net affect on mean bright- 
ness is negligible. Additionally, if our assumptions about the 
dependence of DDT density on metallicity are incorrect and 
further studies indicate the opposite trend, the metallicity ef- 
fect on DDT density could negate the effect due to lepton 
number conservation. 

This study stresses the importance of multi-dimensional ef- 
fects, but has provided evidence of the limitations of 2D sim- 
ulations of supernovae. In the immediate future, we plan 
to perform 3D simulations for better realism. In reality, we 
might expect a two-parameter family in which the morphol- 
ogy of the dominant Rayleigh-Taylor plume is independent of 
the rate of expansion; however, for our 2D simulations, these 
two effects appear to be correlated by the choice of cylindri- 
cal geometry. A plume developing along the symmetry axis is 
like 3D RT compared to one near the equator, which is more 
like 2D RT, and the latter is weaker (slower to develop). In 
addition to developing faster, plumes near the symmetry axis 
represent less volume and, therefore, expand the star less. A 
fast rising plume combined with a slow rate of expansion in- 
dicates there will be a shallow dependence on DDT density 
(since the plume will move through the various DDT densi- 
ties faster and with less expansion). This result may also de- 
pend on our choice that DDT conditions are met at the tops 
of rising plumes. For these reasons, 3D studies are necessary 
to ascertain whether there is a physically motivated correla- 



tion between the dominant plume morphology and the rate of 
expansion and whether this correlation depends on the choice 
of location of the DDT. In transitioning to 3D simulations, we 
expect the growth of many more plumes; however, the rate of 
expansion will likely not exceed that found in 2D. Because 
RT develops faster in 3D, the conditions for DDT will be met 
sooner leading to less overall expansion of the star and more 
mass at high density at the first DDT time. Therefore, we 
expect 3D results similar to the higher yielding 2D realiza- 
tions with a shallower dependence on DDT density. In any 
case, this study indicates that it is possible to determine the 
susceptibility of a particular model to a change of DDT den- 
sity (and, hence, a change in the composition which alters the 
microphysics that set the DDT density). Once the plume mor- 
phology has been established several tenths of seconds into 
the simulated explosion, it should be possible to estimate the 
total ^''Ni yield to fairly high accuracy given DDT conditions. 
However, the treatment of turbulent flame properties may be 
important. 
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